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Abstract. In recent years, there are increasing evidence for a thermal emission component that accompanies the overall non- 
thermal spectra of the prompt emission phase in GRBs. Both the temperature and flux of the thermal emission show a well 
defined temporal behaviour, a broken power law in time. The temperature is nearly constant during the first few seconds, 
afterwards it decays with power law index a ~ 0.7. The thermal flux also decays at late times as a power law with index 
/3 ~ 2.1. This behaviour is very ubiquitous, and was observed in a sample currently containing 32 BATSE bursts. These 
results are naturally explained by considering emission from the photosphere. The photosphere of a relativistically expanding 
plasma wind strongly depends on the angle to the line of sight, 6. As a result, thermal emission can be seen after tens of 
seconds. By introducing probability density function P(r, 9) of a thermal photon to escape the plasma at radius r and angle 
6, the late time behaviour of the flux can be reproduced analytically. During the propagation below the photosphere, thermal 
photons lose energy as a result of the slight misalignment of the scattering electrons velocity vectors, which leads to photon 
comoving energy decay E r (r) °c r~ 2 / 3 . This in turn can explain the decay of the temperature observed at late times. Finally, I 
show that understanding the thermal emission is essential in understanding the high energy, non-thermal spectra. Moreover, 
thermal emission can be used to directly measure the Lorentz factor of the flow and the initial jet radius. 
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INTRODUCTION 

It is widely believed that the prompt emission from GRBs arise from the prompt dissipation of a substantial fraction 
of the bulk kinetic energy of a relativistic outflow, originating from a central compact object. The dissipated energy 
is converted to energetic electrons which produce high energy photons by synchrotron and synchrotron self Compton 
(SSC) scattering. This interpretation was found to be consistent with a large number of GRB observations 0101, which 
generally show a broken power law spectrum at the keV — MeV energy range (which has became known as the "Band 
function" HUSH). 

In spite of its many successes, in recent years there are increasing evidence for low energy spectral slopes that are too 
steep to account for by the optically thin synchrotron - SSC model 10, H,^. Motivated by these findings, an additional 
thermal (blackbody) component was suggested that may contribute to the observed spectrum iflHflllTl^flsll . Indeed, 
from a theoretical point of view, such a component is inevitable: the optical depth near the base of the flow is enormous, 
lye > 10 15 (for a review, see, e.g., [ 14]), thus photons emitted by the inner engine that produces the burst, or by any 
dissipation mechanism that occurs deep enough in the flow, necessarily thermalize before decoupling the plasma 
at the photosphere. While in principle these photons are the first to reach the observer (if emitted on the line of 
sight), in practice, due to the Lorentz contraction the observed time difference between thermal photons originating 
from the photosphere and non-thermal photons originating from dissipation above the photosphere can be shorter 
than millisecond, and thus not be resolved. It should be stressed here, that due to the multiple dissipation processes 
episodes expected during the prompt emission (e.g., by internal shock waves), the existence of thermal photons do not 
contradict emission of non-thermal photons, but adds to it. 

The interpretation of the prompt emission as being composed of thermal emission component in addition to the 
non-thermal one, was put forward by Ryde [15]. In this work, analysis of the time-resolved spectra of nine bright, 
long GRBs which were characterized by hard low energy spectral slopes, showed that a dominant thermal component 
could be used to explain the observed spectra. It was found in this work that the temperature of the thermal component 
is approximately constant and equals a canonical value of T° b « 100 keV for a few seconds, afterwards it decays as a 
power law in time with power law index a ~ 0.6 — 1.1. Ryde suggested later on fl6Tl that a thermal component could 
in fact exist in many bursts in which it does not necessarily dominate over the non-thermal component. 

The suggestion made by Ryde, combined with the theoretical arguments mentioned above, had motivated further 
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FIGURE 1. Left: an example of spectral modeling using a thermal + a single power law to model the non-thermal spectrum. 
Example is shown for BATSE trigger 1974 (GRB921003). Middle and right: temporal evolution of the temperature (middle) and 
flux (right) of the thermal component. The temperature is nearly constant for ~ 1.5 s, afterwards it decreases as a power law in 
time. The flux also shows a broken power law temporal behaviour, with a break time which is within the errors of the break time in 
the temperature evolution. Results are shown for BATSE trigger 3765 (GRB950818). Very similar results are obtained for a large 
sample of bursts (see text). 



work on the origin of the prompt emission, and the role played by the thermal photons. The results of this work are 
published in a series of papers il7lfl8lfl9ll . in which various aspects of the problem are examined. Here, I summarize 
the main results found so far. I first present the key results of the work by Ryde & Pe'er lfl7ll . which provide, for the 
first time, a full analysis (flux and temperature) of the temporal behaviour of the thermal emission component observed 
in a large (32) sample of GRBs. I then present a theoretical interpretation, based on the analysis done in 11811 . Finally, 
I argue that thermal emission, in addition to its contribution (via Compton scattering) to the high-energy, non-thermal 
spectrum, can be used to deduce important physical parameters of the flow, such as the Lorentz factor, F and the initial 
jet radius. 



EVIDENCE FOR THERMAL EMISSION: REPETITIVE BEHAVIOUR 

We use the method developed by Ryde lfl5lfl6ll . to model the prompt emission spectra in the BATSE detector range 
(20keV — 2MeV) as being composed of a thermal component and a single power law. While being incomplete in 
natur^E this model, which contains four free parameter^ (similar to the number of free parameters in the "Band" 
model) provides good statistical fits (x 2 /d.o.f. ~ 1) to the spectra in the BATSE range. An example of this fit is 
presented in figure Q] (left). 

Key results were found when we used this method to model time resolved spectra. In a sub-sample of w 300 long 
BATSE bursts, the spectral and temporal coverages are good enough to enable splitting of the lightcurves into separate 
time bins and model the spectra in each time bin (typically of duration ~ 1/2 s) separately. The first important result 
is that so fafl we were able to identify a thermal emission component, and to model the time resolved spectra in the 
way described above in all the cases studied (currently, our sample contains 32 bursts, in which the thermal component 
does not necessarily dominate the spectra^). By doing so, we retrieved the results found by Rvde lfl5ll (using a smaller 
sample) about the well-defined temporal behaviour of the temperature of the thermal component (see fig. [T] middle): 
the temperature is typically found to be nearly constant for a few seconds, afterwards it decay as a power law in time 
with power law index a ~ 0.4 — 0.9. 

We continued further to analyse the temporal behaviour of the flux of the thermal component (see fig. Q] right). 



A single power law cannot be used to model the data at energies far above or far below the limited BATSE range; nonetheless, theoretically a 
single power law for the non-thermal emission can be justified over a limited energy range. See discussion below. 

2 Temperature and flux of the thermal component, power law index and normalization of the single power law component. 

3 This research is still on-going. 

4 Typically, minimum thermal flux of ~ 10 — 20% of the total flux is needed to be able to identify the thermal component. 
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FIGURE 2. Histograms of the power law indices of the late time decay of the temperature (left) and the thermal flux (right) in a 
sample of 32 bursts. In all the cases, not only we were able to identify a thermal component, but the temporal behaviour of both the 
temperature and flux show a broken power law behaviour, with very similar indices for the different bursts at late times. 



Here, too, we found a well defined temporal behaviour: the thermal flux slightly rises for a few seconds, afterwards 
it decays as a power law in time with power law index j3 ~ 2.1. In all the cases studied, the break time in the flux is 
within the errors of the break time in the temperature. The most important result is the repetitive behaviour of both 
the temperature and the flux of the thermal component: a similar temporal behaviour (broken power law of both the 
temperature and thermal flux) was found in all the cases studied so far, with very similar power law indices. Histograms 
of the late time power law indices for the sample of 32 bursts are presented in figure [2] 



THEORETICAL MODEL AND EXPECTATIONS 

Our basic assumption in an attempt to understand these results is that thermal photons originate either from the inner 
engine that produces the relativistic outflow, or from an unspecified dissipation process that occurs deep enough in 
the flow, so that the photons thermalize before escaping the plasma at the photosphere0. The angular dependence of 
the photospheric radius in a relativistically expanding plasma, characterized by constant Lorentz factor F and constant 
mass ejection rate was first studied by lEofl . Recently, we showed [18] that it can be formulated in a surprisingly simple 
form, 
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This function is plotted in figure|3](left) for two values of F. Here, R c / is a constant which depends on the mass ejection 
rate, 9 is the angle to the line of sight, j3 is the plasma expansion velocity and the last equality holds for 9 -C 1, r> 1. 
The strong 0-dependence implies that for characteristic GRB luminosity L = 10 52 Lergs _1 and F = 10 2 r2, thermal 
photons escaping the photosphere from high angles to the line of sight 9 w 0. 1 i (estimated GRB jet opening angle; 
see, e.g.. lf2lll ). are delayed with respect to photons originating on the line of sight by rs 30Lr2 1 0ij s. 

The definition of a photosphere as a surface in space from which the optical depth to scattering equals unity, is 
however, incomplete: photons have a finite probability of being scattered at every point in space in which electrons 
exist. Therefore, in order to fully quantify the last scattering event positions, one needs to use probability density 
function P(r, 9). Using the simplified assumptions that the last scattering event radius is independent on the scattering 
angle, and that in the (local) comoving frame the scattering is isotropic, we showed ( 11811 ) that the probability density 
function can be written as 

p ^ = {j) 2r 2 /3[i-cos(e)] 2 ' (2) 



5 A possible alternative model is emission of thermal radiation by dissipation above the photosphere. 
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FIGURE 3. Left: Angular dependence of the photospheric radius rp/,(0;r) for T = 5, 100 (see eq.QJ. Right: result of a Monte- 
Carlo simulation shows the positions of the last scattering events of photons originating from below the photosphere. The green 
line shows rp/,(r = 100), which by definition is the surface in space from which the optical depth iy e = 1. Contour lines are added 
to show the distribution of the last scattering events positions. Clearly, last scattering events positions are located (with different 
probabilities) in the entire space, while r p i, as defined here only gives a first approximation. 



where ro = r p /,(0 = 0) = Rj / (2kF 2 ). In order to validate the approximations used, as well as the assumptions of the 
diffusion model presented below, we carried a Monte Carlo simulation that traces the photons from deep inside the 
flow until the final scattering event. The results of this simulation are presented in figure [3] (right). 

The probability density function enables an analytic calculation of the late time decay of the flux and energy of the 
thermal photons originating from below the photosphere. In the calculation, we use the diffusion approximation, in 
which all the photons are injected into the flow at the center of the expansion, at time t = (i.e., 5-function injected 
in space and time). The photons are coupled to the flow until the last scattering event takes place. Therefore, before 
decoupling, the velocity component of the photons in the direction of the flow is w j3c. The observed time delay of 
a photon whos last scattering event is at (r, 0) with respect to a "trigger" photon that was emitted at the center of the 
flow at t = and was not scattered at all, is thus t ob — At" b = (r/j3c) [1 — j3 cos(0)]. The observed flux is calculated by 
integrating the probability of a photon to be scattered over the entire space, while maintaining the correct arrival time, 
F" b {t" b ) oc Jdrfd9P(r, 9)8 (t ob = (r/j8c)[l - J3cos(0)]). At late times, this gives F ob {t" h ) °c t ob ~ 2 (see fig.g] left). 

The observed energy of a photon is blue shifted by the Doppler factor &{9) = [r(l — j3 cos(0)]~' with respect to 
its (local) comoving energy, which itself depends on the photon propagation radius within the flow, e' = e'(r), via 
two effects: the first is adiabatic energy losses of the scattering electrons. The second is energy loss of the photon 
due to the slight misalignment of the scattering electrons velocity vectors. Thus, even if there is no energy exchange 
between an electron and a photon in a single scattering event (i.e., Thompson scattering), the next scatterer's velocity 
vector is not parallel to the first ones, hence the photons' energy in the frame of the next scatterer is slightly lower. 
We showed in 111811 , that deep inside the flow this effect leads to photon (local) comoving energy loss, e'(r) oc r - 2 P/ 3 ^ 
which relaxes as the photon approaches the photospheric radius ro to e'(r;r > ro) oc r°0 Using again the probability 
density function defined in equation [2] the temporal evolution of the observed energy of photons originating from 
below the photosphere is obtained by an integration over the entire space in a similar way to the calculation of the 
flux, e" b {t" b ) oc Jdrfd9P(r,9)e'(r)3i(9)8 (t ob = (r/j3c)[l - j3cos(0)]). At late times, this gives e ob {t" b ) oc t" b ~ a , 
with a w 1/2 — 2/3. The analytical results, together with the more accurate numerical results are plotted in figure [4] 
(middle). 

The theoretical model thus gives well-defined predictions for the late time thermal flux and temperature decay: 
provided that the inner engine decays fast enough, the observed flux of thermal photons is expected to decay as t ~ 2 , 
and the observed temperaturfl as t~ a with a ~ 1/2 — 2/3. The characteristic time scale of the decay is predicted to 



6 This effect is very similar to energy loss by adiabatic expansion of the photons. Note though that the conditions here are somewhat different than 
that of classical adiabatic expansion, since the photon propagation volume is in principle not limited. 

7 In the analysis presented we treated single photons, while a temperature is defined for Plank distribution of photons. Although the observed 
spectrum deviates from Plank spectrum, being a convolution of Planck spectra it is not expected to deviate much from it. 




FIGURE 4. Left: observed thermal flux as a function of the observed time, presented in units of normalized time, f# = 
ro(l — j8)/c. The blue (solid) line is the numerical simulation result, and the red (dash) line is the analytical approximation. 

Clearly, the analytical approximation is good at late times, where the thermal flux decays as F" b (t ob ) <>= t" b . Middle: same 
for the observed temperature. At late times, the simulation results and the analytical model predict observed temperature decay 
T ob {t ob ) <x t ob a , with a ~ 1/2 — 2/3. Right: examples of time averaged spectra obtained for different values of the optical depth 
ly e at the dissipation radius, under the assumption that thermal component exists. Ty e = 1 represents dissipation at the photospheric 
radius, r$. The electrons energy distribution is modified by multiple Compton scattering with the thermal photons, and as a result 
the spectrum deviates from a broken power law typical for optically thin synchrotron- SSC emission. See 1 23, 24] for details. 



be tens of seconds. 



IMPLICATIONS: COMPTON SCATTERING AND PROPERTIES OF THE FLOW 

The existence of thermal emission component is potentially crucial in understanding not only the spectrum in the 
BATSE range, but the very high energy (> MeV) spectra as well. Thermal photons can serve as se ed p hotons for 
Compton scattering by energetic electrons produced by dissipation processes in the flow lfl3l I22I l23i l24ll . Since the 
nature of the dissipation processes (e.g., internal shock waves or magnetic reconnection) is yet unclear, it may occur 
at a variety of radii, including near or below the photosphere. In this case, energy exchange via both inverse and 
direct Compton scattering with the thermal photons may significantly modify the electrons energy distribution, and 
as a consequence a variety of non-thermal spectra, which cannot be described by a simple broken power law may be 
obtained I23l 12411 . The effect on the high energy spectra may be significant even if the dissipation occurs at radii which 
are 1-2 orders of magnitude above the photosphere. Examples of possible spectras are plotted in figure [4] (right). 

In addition to their role as seed photons for Compton scattering, thermal photons can be used to directly probe 
the properties of the flow. The theoretical model presented above is able to explain the late time temporal decay of 
both the temperature and the flux of the thermal emission. If this explanation is correct, it implies that the thermal 
photons observed at early times (before the temporal break) are emitted on the line of sight. The dimensionless ratio of 

the thermal flux and temperature M = (F ob / ctT * 4 ) 1 / 2 , where a is Stefan's constant, is thus proportional to ro/Fdi, 
where di is the luminosity distance, and the Lorentz factor T originates from relativistic aberration. For constant flow 
velocity, ro = r p /,(0 = 0) °c L/r 3 (e.g., [25']), where the luminosity L can be determined once the flux and the redshift 
are known. Thus, both the Lorentz factor T (L/ Mdt) 1 ^ and the photospheric radius ro can be directly determined 
for bursts with known redshift and identifiable thermal component. 

In principle, ro is the innermost radius from which information can reach the observer. However, the fireball model 
predicts the dynamics of the plasma below the photosphere using energy and entropy conservation (e.g., j25lP . Using 
these assumptions, we showed ( lfl9lo that the base of the je^is °< d^, and thus can be determined. Using this method 
for GRB970828 at redshift z = 0.96 we found T = 305 ±28 and r base = (2.9 ± 1.8) x 10 s cm. These results are 
consistent with earlier estimates, based on light crossing time arguments and early afterglow emission measurements. 
Former methods, though, can provide either lower limit or values estimates good to an order of magnitude, while the 



8 Defined here as the radius at which T = 1 ; rb ase may also be identified with the sonic point. 



statistical error in the method presented here on the value of T is ~ 

SUMMARY 

In this work we examine various aspects of thermal emission from GRBs. Ryde & Pe'er [17] use the method suggested 
by Ryde llaJlal to analyze the prompt emission spectra in an alternative way to the commonly used broken power law 
(the "Band" model). By doing so, we introduce a new physical meaning to the spectrum, as being composed of thermal 
+ non- thermal emission. We find a repetitive behavior of both the temperature and the thermal flux, in a sample of 
32 bursts. The late time decay indices of the temperature and flux are consistent with the predictions of the theoretical 
model developed by Pe'er 01811 . based on the idea of extended photospheric emission from higher angles and higher 
radi{3- We showed that thermal emission must be considered in order to correctly interpret emission at higher energies 
(at the GLAST energy band). Moreover, we showed ([ 19]) that thermal emission can be used to determine important 
parameters of the GRB outflow, such as the Lorentz factor and the radius at the base of the jet. 

We find the repetitive behaviour and the agreement between the theory and observations very encouraging. We 
continue our work on this project, as we believe that it could provide new understanding of the mechanism and physics 
of the prompt emission and of GRB progenitors. 
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9 Systematic uncertainty due to the unknown non-thermal flux also exists; this though can in principle be removed using late time measurements. 

10 This idea has some similarities to the "high latitude emission" models used to model GRB afterglow emission. However, here we treat emission 
from optically thick rather than optically thin plasmas. 



